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ABSTRACT 

We investigate two successive flux rope (FRl and FR2) eruptions resulting in two coronal mass 
ejections (CMEs) on 2012 January 23. Both FRs appeared as an EUV channel structure in the images 
of high temperature passbands of the Atmospheric Imaging Assembly prior to the CME eruption. 
Through fltting their height evolution with a function consisting of linear and exponential components, 
we determine the onset time of the FR impulsive acceleration with high temporal accuracy for the flrst 
time. Using this onset time, we divide the evolution of the FRs in the low corona into two phases: a 
slow rise phase and an impulsive acceleration phase. In the slow rise phase of the FRl, the appearance 
of sporadic EUV and UV brightening and the strong shearing along the polarity inverse line indicates 
that the quasi-separatrix-layer reconnection likely initiates the slow rise. On the other hand for the 
FR2, we mainly contribute its slow rise to the FRl eruption, which partially opened the overlying 
fleld and thus decreased the magnetic restriction. At the onset of the impulsive acceleration phase, 
the FRl (FR2) reaches the critical height of 84.4±11.2 Mm (86.2±13.0 Mm) where the decline of the 
overlying fleld with height is fast enough to trigger the torus instability. After a very short interval (~2 
minutes), the flare emission began to enhance. These results reveal the compound activity involving 
multiple magnetic FRs and further suggest that the ideal torus instability probably plays the essential 
role of initiating the impulsive acceleration of CMEs. 

Subject headings: Sun: corona — Sun: coronal mass ejections (CMEs) — Sun: flares — Sun: magnetic 
topology 
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1. INTRODUCTION 

Coronal mass ejections (CMEs) are the most spectac- 
ular eruptive phenomena in our solar system. They are 
able to release a large quantity of plasma and magnetic 
flux into the heliosphere and sever ely disturb t he space 
environment around the Earth (e.g., Gosling|199 3). Over 
the last 20 years, although the solar community has made 
a considerable progress in many aspects of understand- 
ing CMEs, the i mportant issue of how CMEs are initi ated 



initiating the explosive er uption of CMEs. Differing from 



remains elusive (Chen 2011 Schmieder et al. 2012). In 



this Letter, we report the compound eruption activity 
involving two inter-connected flux ropes (FRs; magnetic 
fleld lines wound around each other) in the same active 
region, and then elucidate their initiation mechanism. 

In terms of whether involving magnetic reconnec- 
tion, existing initiation models can be divided into two 
groups; one group are reconnecti on-based m odels includ- 
ing tether-cutting recon nection ([Moore et al. 2001) and 
breakout reconnection ([An tiochos et al.||1999| ^Chen & 
OOOl iKarpen et al. 2012), and the other 



2000| [Karpen 
Dasea ideal magnetohydrodynamics (M HD) mod- 



group 



Shibata 
areTTT 

els, involving c atastrophic loss-o f-equilibrium (Forbes 
ttsenberg 199l]_ Isenberg et al.| p^993), ki nk instability 
([Torok et al.|2QQ4| ), and torus instability ( |Kliem fc Torok 
2006p . In the tether-cutting (breakout) model, key mech- 
anism solely concerns the magnetic reconnection in the 
CME core (overlying) fleld region that increases upward 
magnetic pressure (reduce the overlying tension), thus 
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the reconnection mo dels, Forbes & Isenberg (1991) and 



Isenberg et al.|(^1993) showed that the 1^'K will lose equi- 



librium in an ideal MHD process if the photospheric 
sources of the overlying fleld converge toward each other. 
In addition to this catastrophic loss of equilibrium, kink 
or torus instability is also capable of initiating the CME 
explosive eruption. The ideal kink instability develops if 
the averag e twist of the FR is greater than a threshold 
(e.g., 3.57r; Torok et al.|2004 ). The torus instability takes 
place if the restoring force caused by overlying fleld de- 
creases faster than the outward dir ected Lorenz self-force 
as the flux rope expands outward ( Kliem fc TQrQk|[2006 
Olmedo fc Zhang_2010 



/ith the development of these theoretical models, val- 
idating and distinguishing them observationally becomes 
a matter of great necessity. We here investigate in detail 
the initiation of a compound CME activity originating 
from two successive FR eruptions observed by the At- 
mospheric Imaging Assembly (AIA; jLemen et al. 2012|) 
on board Solar Dynamics Observatory [SDO). We find 
that the slow rise of the flrst CME is most likely due 
to the quasi-separatrix-layer (QSL) reconnection in the 
low corona; while the slow rise of the second one results 
from the partial opening of the overlying fleld by the flrst 
CME. However, for their initiation of the impulsive ac- 
celeration, we contribute the mechanism to be the ideal 
torus instability. Data and results are presented in Sec- 
tion 2 and 3, respectively. Summary and discussions are 
given in Section 4. 



2. INSTRUMENTS 

The data we used are mainly from the AI A, as weh as 
the H ehoseismic and Magnetic Imager (HMI; Schou et al. 
2Q12D on broad SDO. AIA provides the EUV images of 
the solar corona with the temporal cadence of 12 sec- 
onds, the pixel size of 0.6'^ and the field of view (FOV) 
of 1.3i?0. The HMI provides the vector magnetic field 
of the solar photosphere with almost the same pixel size 
and FOV as AIA but the temporal cadence of 12 min- 
utes. Th e Large Angle and Sp ectrometric Coronagraph 



Brueckner et al.|1995 ) on board SOHO and the 
1 Connection Coronal and Heliospheric Investi- 



(LASCO; 

Sun-Eart 

gation (SECCHI; "Howard et al."2008) on board STEREO 
observe the CMEs in white light. GOES X-ray data re- 
veal the SXR (SXR) 1-8 A fiux of CME-associated fiares. 

3. RESULTS 

3.1. Kinematics of CMEs 

On 2012 January 23, there are two CMEs (CMEl 
and CME2) that appeared successively in the FOVs of 
LASC0/C2 and SECCHI/CORl, finally manifesting as 
a compound CME in the FOV of LASC0/C3 (als o see; 



|Joshi et al.||2Q13[|Shen et al ||2Q13[ [Liu et alj2Q13[ ). The 

selected snapshots are displayed in Figure [l| Here, we 



take advantage of th e graduated cylindrical s hell (GCS) 
model developed by Thernisien et al. (2006) to recon- 
struct the three-dimensional (3D) morphology of the 
CMEs and obtain the true heights of the CMEs. The 
forward fitting model includes three positional parame- 
ters and three dimensional parameters: the Carrington 
longitude ((/)) and latitude (^), the tilt angle (7) with re- 
spect to the equator, the height {h) and aspect ratio {n) 
of the FR, and the half-angle (a) between the two FR 
legs. For the CMEl at -03:30 UT, the fitting parameters 
are ^=219°, 6>=25°, 7=-49°, h = 4.1 Rq, /^=0.4, and 
a=32°; for the CME2 at -04:00 UT, ^=210°, 6>=32°, 
7=-71°, h = 4.2 Rq, a^=0.4, and a=40°. The corre- 
sponding wireframe rendering are shown in Figure [l] 

The FR of CMEs have been repor ted to appear as an 
EUV hot channel structure (H CS) jZha ng et al.||2012[ 
Cheng et al. (20121 [20131 | Li fc Z hang 2013) and/or bulplple 



Cheng et al. '201l[ Su et al.[,2012| |Patsourakos et al. 



013[) in the AIA images, bhi the two CMEs studied 



we inspect the AIA images and find that the FRl only 
appeared in 131 A and 94 A passbands (Figure ^b)) 
but not in any other ones, indicating that it should be 
hot(>7MK;[0'Dwyer et al1[2Q 10, Reeves fc Golub|[2QTT 



Cheng et al. 201l| ). On the other hand, the FR2 is visible 
in all AIA passbands, although most noticeable in 131 
A and 94 A images (Figure M^e)). The mixture of hot 
and cold plasmas within the FI12 probably is attributed 
to th e presence of filament material (e.g. 



2012[ ) . More detailed evolution of the two HCSs please 
refers to the attached movies. In order to clearly display 



the rising motion of the HCSs, as well as the expansion of 
the overlying field, we make two slices (slicel and slice2 in 
Figure |2|b) and (e)) alongthe HCS rising directions and 
one slice (slice3 in Figure [2|h)) along the overlying field 
expanding direction. The time evolution of the HCSs 
and the overlying field along these slices are shown in 
the slice-time plots in Figure [sF a) -(c). Through these 
time-stacking plots, we measure the height of the two 
HCSs with time (diamonds in Figure |3ja) and (b)). The 



height-time data are plotted in Figure[3[d). From Figure 
S[g)-(i), as well as the stack plot of the slice3 (as shown 
Toy two inclined lines in Figure IsFc)), we notice the slow 
expansion of the overlying fieloafter the FRl eruption, 
which most probably leads to the initial rise motion of 
the FR2. 

With the height-time data, we further calculate the 
velocity and acceleration. In order to reduce the uncer- 
tainty in the height measurement, a cubic spline smooth- 
ing method is used t o smooth the height (also see, Pat-l 
sourakos et al.| 2010 Cheng et al.^2013} . Then a piece- 
wise numerical de rivative method is applied to calculate 
the v elocity (e.g., 'Zha ng et aL]|2001[ |2Q04| [Cheng et al^ 
2010[ ). The deduced velocity-time profile is displayed in 
Figure |3|e). The uncertainty in the velocity is mainly 
from the uncertainty in the height measurement (to be 4 
pixels; 1700 km for AIA observations and 44,000 km for 
GCS fitting). Similarly, we derive the CME acceleration 
and resulting uncertainty, as shown in Figure [sFf ) . Note 
that all heights refer to the top of the CME FRs from 
the solar surface. 

From Figure IsFe) and (f), one can find that the kine- 
matical evolution of the CMEs are tightly associated with 
the emission variation of associated fiares. The velocity 
evolutions of the CMEl and CME2 are very consistent 
in time with the SXR 1-8 A fiux profiles associated with 
fiarel and fiare2. The acceleration of the CMEl and 
CME2 grows in step with the derivative of the SXR pro- 
file. The close temporal correlation between the CME 
kinematics and the fiare emission implies that CMEs and 
associated fiares are not two independent eruption phe- 
nomena but only t wo distinct manifestations of the same 
eruption process (|Lin fc Forbes 2000| Priest fc Forbes 



20021 [Zhang fc Dere||2006| [Temmer et al.||2010| ) 



3.2. Onset of Impulsive Acceleration of CMEs 

In order to investigate the initiation of the CMEs, we 
first determine the onset of the impulsive acceleration. 
Assuming that the height evolution of the CME FR in 
the low corona follows a function h{t) = coe*^^~^°^/'^ + 
ci(t— to)+C2, where h{t) is height, t is time, r, to, cq, Ci, C2 
are five free coefficients. The model includes two compo- 
nents: the linear and the exponential, which correspond 
to the slow rise phase characterized by a constant veloc- 
ity and the initial impulsive acceleration phase character- 
ized by an exponentially increase of velocity, respectively. 
The model is reasonable due to the fact that the velocity 
of the CMEs increases rapidly once the impulsive ac- 
celer ation is triggered whether by the flare reconnection 




Antiochos et al. 1999 



[Moore et al.| [2001j jkarpen 
or by the torus i nstability (e.g., [Torok & 
Olmedo fc Zhang[[2010|. The fitting to~tIie 



Cheng et al. oat a is achieved by the "mpfit.pro" routine. With the 

fitting parameters, the onset of the CMEs is defined as 
a time where the exponential velocity is equal to the lin- 
ear velocity (tonset='?"ln(cir/co) +to). Accordingly, the 
height at the onset time corresponds to the critical height 
h{tonset) of the eruption. 

We further use 100 Monte Carlo (MC) simulations to 
estimate the uncertainties of our results. For each MC 
realization, the measured heights are first perturbed ran- 
domly by an amount of 5 within a sigma equal to the 
error of the height, and then re-fitting the model to the 
heights. The final onset time and onset height are the 



averages of 100 onset times and onset heights derived 
by 100 MC reahzations. The corresponding uncertain- 
ties are triple the standard deviations {3a) of the MC 
fitting. The fitting results are shown in Figure [4j For 
the CMEl, we determine the onset at 02:02 UT with the 
error of 2 minutes; the onset height is 84.4±4.2 Mm. For 
the CME2, the onset is 03:34 UT with the error of 1 
minute; the onset height is 86.2±6.0 Mm. Moreover, the 
uncertainty of the reference point of measuring height is 
estimated to be ~7.0 Mm (10^^). Therefore, the final un- 
certainties of the onset heights for the CMEl and CME2 
are '^11.2 Mm and 13.0 Mm, respectively. 

From Figure [4Fc) and (f), one can find that the on- 
set of the CMEs is earlier than that of associated flares 
by a few minutes. Here we define the onset of flares as 
a time where the derivative of the SXR flux becomes 
positive and begins to increase successively. For CMEl 
and CME2, the leading times, i.e., the onset of CME 
impulsive acceleration relative to the onset of the flare 
impulsive phase, are about 2 minutes. Actually, for the 
flare2 the onset time in GOES record is 03:38 UT; the 
leading time would be of ~4 minutes if this onset time is 
adopted. The results imply that the impulsive accelera- 
tion onset of the CME FRs is most probably caused by 
ideal MHD instability rather than by the flare reconnec- 
tion. In addition, we note that the impulsive acceleration 
of the FRs should occur earlier than that we determined 
if the projection effect for the heights in the AIA FOV 
is taken into account. The result should be strengthened 
because it is even earlier for the FRs to ascend to the 
critical height of the torus instability. 

3.3. Magnetic Properties of CMEs 

Previous theoretical works have proved that the oc- 
currence of torus instability of a FR depends on the 
decay index of the backgro und magn etic field B with 
height h (n=-d\iiB/d\nh;\K\iem fc T orok 2006; Isen- 
berg fc Forbesl2007l[Olmedo^ Zhang|2010nDemoulin fc, 



i. 



Aulanier]2010p . The decay index is computed from the 
potential field model based on the vertical component 
of vector magnetograms provided by the HMI (Figure 
5|a)-(c)). The distributions of the average decay index 
with height over three regions (the main polarity inverse 
line (PIL), two rectangle regions near the main PIL; Fig- 
ure pjc)) are shown in Figure |5|d). One can find that 
at tne onset of the CME impulsive eruption, the FRl 
(FR2) reached the height of 84.4±11.2 Mm (86.2±13.0 
Mm) where the decay index is 1.7±0.1, which is slightly 
larger than the nomin al critical value^l.5 of the torus 
instability occurrence (Torok & Kliem'2005'). This result 
supports the theory of the torus instability as the trigger 
of the impulsive acceleration of CMEs. Note that here 
the decay index is calculated from the horizontal com- 
ponent of 3D magnetic field because the vertical compo- 
nent does not have a role in constraining the FR. Also 
note that, in torus instability models, the field induced 
by the current inside the FR and the background field 
constraining the FR constitute the total magnetic field, 
which usually can be modeled by the nonlinear force- fre e 



field (NLFFF) (e.g., ^Guo et al. 2010 Sun et al. 2012) 



Considering the fact that it is difficult to separate the 
background field from the FR in the NLFFF model, the 
potential field thus is used to be an approximation of the 
background field (also see. Fan fc Gibson||2007 Aulanier 



iriI]|2QTQl [Cheng etal.|[2QT3 ). 

Based on the HMI vector data at 00:00 UT on Jan- 
uary 23 (Figure |5jc)), we extrapolate the 3D NLFFF 
using the optimization algorithm in cartes ian coordinate 
(Wheatland et al."2 000[[Wie gelman n|2004p . The selected 
held lines are displayed mTlgure Isfe). One can find that 
the extrapolation suggests there exist two FR structures 
(marked as FRa and FRb), both of which have an associ- 
ated filament visible in the AIA 304 A passband (Figure 
[5lf)). The FRa finally erupted as the CME2, while the 
FRb always stayed there during the whole eruption pro- 
cess. However for the FRl discussed later, we are not 
able to reconstruct a corresponding FR even resorting to 
the optimization algorithm in spherical geometry with a 
larger FOV (e.g., Guo et al.|2012[ ). This may not be sur- 
prising, since the FRl as observed in AIA is a larger and 
higher structure, which can not be well modeled with the 
current extrapolation technique. 

4. SUMMARY AND DISCUSSIONS 

In this Letter, we report the initiation process of a 
compound CME activity consisting of two successfully 
erupted FRs. We find that the kinematics of the FRs in 
the low corona have two phases: a slow rise phase and 
an impulsive acceleration phase, which can be fitted very 
well by a model consisting of a linear component and an 
exponential component. 

In the slow rise phase of the FRl, we inspect the AIA 
images in all EUV and UV passbands and find that some 
brightenings spread sporadically along the PIL under 
the FRl and at its two footpoints as well. It indicates 
that magnetic reconnection probably take place inside or 
around the FRl. The reconnection is most likel y to be at 



the QSL surrounding the FR in the low corona rtAulanier 



ik 



jCt al.j,2010p . The QSL reconnection is too weak to prO' 
duce nonthermal particles. It is different from the recon- 
nection process occurring during the flare, in terms of 
both intensity and location. With this QSL reconnection 
taking place, the FRl is allowed to grow and rise slowly. 
However for the FR2, the driving mechanism of the slow 
rise phase might be different from that of FRl. Due to 
the eruption of the FRl, the overlying field is opened par- 
tially, leading to the expansion of the ambient magnetic 
field. The expansion decreases the downw ard magnetic 
tension, causing the slow rise of the FR2 (|Torok et al. 
[2QTT| [Lynch fc Edmon ds"on|[2Q1 3| 



Through 3D NLFFF extrapolation, it is indicated that 
the FR2 exists prior to the slow rise phase for a long 
time. While for the FRl, we do not know whether it 
exists or not prior to the slow rise. It may be formed 
locally through the QSL reconnection in the slow rise 
phase (e.g., I Aulanier et al.|2010[ Liu et al.|2010 ) or there 
has existed a nascent bR, and the QSL reconnection only 
plays a role in enhancing the poloidal flux of the FRl. 

In the end of the slow rise phase, the FRl (FR2) 
reaches the height of 84.4±11.2 Mm (86.2±13 Mm), 
where the decrease of the background magnetic field with 
height is fast enough such that the torus instability takes 
place, thus triggering the impulsive acceleration. It is 
worth mentioning that we do not find evidence of the 
kink instability from the observations. 

As the impulsive acceleration commences, the FR 
quickly stretches the anchored overlying field upward, 
and then the magnetic field underneath starts to recon- 



nect impulsively. Such runaway reconnection is the cause 
of the observed flare rise phase. The time lag between 
the onset of associated flare and that of torus instability 
is only a few minutes, probably shorter. Without high 
cadence observations from AIA, it would be difficult to 
distinguish the relative timing between flux rope impul- 
sive acceleration and flare reconnection onset. 

As the flare starts, the CME acceleration and the flare 
reconnection are coupled together. The flare reconnec- 
tion is able to rapidly convert ambient magnetic field 
into the FR, enhancing the upward Lorentz self- force, 
thus further accelerating the FR. In response to the es- 
caping of the FR, reduced magnetic pressure would drive 
a plasma inflow, which in turn causes more ambient mag- 



netic field reconnecting. Therefore, the runaway recon- 
nection after the flare onset is in a positive feedback pro- 
cess that impulsively accelerates the CME and enhances 
the flare emission. 
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Fig. 1. — White light coronagraph images and GCS FR (blue and green lines) modeling of the CMEl and CME2 on 2012 January 23. 
Blue and green circles denote the location of the solar limb. 
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Fig. 2. (a)-(f) AIA 94 A (-7 MK) images of the CMEl and CME2 HCS. (g)-(i) AIA 193 A (-1 MK) images of the overlying field. In 
panels (b), (e), and (h), overlaid contours show the positive (purple) and negative (black) polarity of magnetic field. The two arrows point 
out the two FRs. The three oblique solid lines indicate the orientations of the slicel, slice2, and slice3. 

(Animations of this figure are available in the online journal.) 
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Fig. 3. — (a)-(c) Stack plots of the intensity along the slicel, slice2, and slice3 as shown in Figure [2] The diamonds in panel (a) and (b) 
indicate the height-time measurements of the CME HCSs, two vertical lines in panel (c) denote the onsets of the FRl and FR2, two inclined 
lines show the expansion of the overlying field, (d)-(f) Temporal evolution of the heights, velocities, and accelerations of the CMEl and 
CME2. The GOES SXR 1-8 A fiux and resulting time derivation are shown in panels (e) and (f), respectively. Note that the acceleration 
of the FRl is tripled to compare with the time derivation of the fiarel. 
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Fig. 4. — (a)-(c) Temporal evolution of the height, velocity, and acceleration of the CMEl HCS in the FOV of AIA. The black solid lines 
show the model fitting to measured height, velocity, and acceleration, respectively. The red solid lines show the GOES SXR 1-8 A fiux and 
resulting time derivation. The vertical blue (horizontal) line points out the onset time (height) of the FRl. The vertical red line indicates 
the onset time of the fiarel. (d)-(f) Same as (a)-(c) but for the CME2 HCS. The vertical green (horizontal) line denotes the onset time 
(height) of the FR2. The vertical red line shows the onset time of the fiare2. 
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Fig. 5. — (a)-(c) HMI vector magnetograms with the 180° ambiguity removed. The background images are the processed vertical 
magnetograms. The green boxes indicate the main AR that the two HCSs origin from. The red curve in panel (c) shows the main PIL. The 
two black boxes point out two regions (Rl and R2) near the main PIL. (d) Distributions of the horizontal magnetic field decay index with 
height over the main PIL, Rl and R2. The two vertical lines display the onset heights of the CMEl and CME2, respectively. The widths 
of the blue and green bars show the uncertainties in the onset heights, (e) Top view of the extrapolated 3D magnetic field configuration. 
The background image is the processed vertical magnetogram overlaid with the contours. The two arrows indicate the FRa and FRb- (f) 
AIA 304 A (^^0.05 MK) image overlaid with the contours of line-of-sight magnetic field, the purple (black) means the positive (negative) 
polarity. The two arrows denote the filament A (Fa) and filament B (Fb). 

(An animation of this figure is available in the online journal.) 



